I 

a: 

o : 
^ ; 

O 

Oh: 



> ■ 
^ ■ 
o ; 
o 
vn ■ 
o ■ 
' 

ON . 

o ; 
•I— I 

^ ; 



X 



Low-energy elementary excitations of a trapped Bose-condensed gas 
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We develop the method of finding analytical solutions of the Bogolyubov-De Gennes equations for the excitations of 
a Bose condensate in the Thomas- Fermi regime in harmonic traps of any asymmetry and introduce a classification 
of eigenstates. In the case of cylindrical symmetry we emphasize the presence of an accidental degeneracy in the 
excitation spectrum at certain values of the projection of orbital angular momentum on the symmetry axis and 
discuss possible consequences of the degeneracy in the context of new signatures of Bose-Einstein condensation. 
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The recent realization of Bose-Einstein condensation 
(BEG) in trapped alkali atom gases 0-0 , followed by 
the second generation of experiments y-]?], has opened 
the possibility of investigating macroscopic quantum phe- 
nomena in these systems. For understanding the macro- 
scopic quantum behavior of a trapped Bose-condensed 
gas especially important is the character of elemen- 
tary excitations of the trapped condensate, which to a 
large extent is predetermined by the interaction between 
atoms. In dilute gases the interaction is primarily bi- 
nary and is characterized by a single parameter, a, the 
s-wave scattering length. This allows one to develop a 
transparent theory which can be tested experimentally. 

At present theoretical investigations of elementary ex- 
citations of trapped Bose condensates include analyti- 
cal solutions for the spectrum of low-energy excitations 
in spherically symmetric harmonic traps in the Thomas- 
Fermi regime |^ and numerical analysis of the eigenfunc- 
tions and eigenenergies of the excitations in the traps of 
spherical and cylindrical symmetry jl^-Q. In the lat- 
ter case the eigenfrequencies of the lowest excitations, as 
those measured in the JILA and MITJ zl experi- 
ments, have also been found analytically |p|, p^ l^Jl^ ■ 

Most interesting are the low-energy excitations, i.e., 
the excitations with energies much smaller than the 
chemical potential (mean field interaction between parti- 
cles), as they are essentially of collective character. Pre- 
vious studies revealed that the eigenfrequencies of con- 
densate oscillations are strongly different from those of a 
collisionless thermal gas , but are rather close to the 

frequencies of a thermal gas in the hydrodynamic regime 
[^,0. In this paper we develop the method of finding 
analytical solutions of the Bogolyubov-De Gennes equa- 
tions for the spectrum and wavefunctions of the conden- 
sate excitations in the Thomas- Fermi regime in harmonic 
traps of any type of asymmetry and introduce a classi- 
fication of eigenstates. We analyse the structure of the 
excitation spectrum in the case of cylindrical symmetry 
and find an accidental degeneracy at certain values of the 
projection of orbital angular momentum on the symme- 
try axis. We address the question of how the accidental 



degeneracy can manifest itself, providing us with a clear 
distinction between the condensate oscillations and the 
oscillations of a classical gas in the hydrodynamic regime. 

We consider a Bose-condensed gas in an external har- 
monic potential V{r) — M^-ujfrf/2 with frequencies 
LUi and assume a pair potential of the atom-atom interac- 
tion of the form U(R) = US(R), where U = Anh^a/M, a 
is the scattering length and M the atom mass. Then the 
grand canonical Hamiltonian of the system is written as 



H=Jdr'^\r] 



-A_A + l/(r)-^ + i;7^t(r)^(r) 



f (r), (1) 



where ^(r) is the field operator of atoms, and /x the 
chemical potential. The field operator can be repre- 
sented as a sum of the above-condensate part and the 
condensate wavefunction = ('!/), which is a c-number: 
= ^' -I- ^fp (see Assuming that the condensate 

density greatly exceeds the density of above-condensate 
particles we omit the terms proportional to and 'F'^ 
in Eq.(|^) and write the Hamiltonian in the form 



H^Hq+ dr< 



r|^'^(r 



2M 



A + Vir)-n 



4|*o(r)|'x 



^r't (r)^'(r) + (i.)^'t (r) + ^-*2^/ (i.)^.'(r)] | ; 



(2) 



Ho^ dr%{r) 



-l-A + V{v)-^, + ^U\Mr)\' 



^o(r). (3) 



The Gross-Pitaevskii equation for ^Fq normalized by the 
condition / |5'o(r)p(ir ~ Nq {Nq is the number of parti- 
cles in the condensate) follows directly from Hq in Eq.(0) 



A + V(r) 

2M ^ ' 



(4) 



Owing to Eq.(^) the part of the Hamiltonian, linear in 
^' (and not included in Eq.(||)), becomes equal to zero. 
The Hamiltonian (^) is bilinear in the operators 
and can be reduced to a diagonal form 



H = Ho + J2 ^Mbu 



(5) 



1 



by using the Bogolyubov transformation generalized to 
an inhomogeneous case: 5''(r) =^^[{uu{'r)K — v*{r)bl] 
Ip^ . Here b^, b\, are creation and annihilation operators 
of elementary excitations. The Hamiltonian takes the 
form (H) if the functions Vi, satisfy the equations 



-2]gA + V^(r)jM, + C/|*o|'(2w.-«.) = (// + ^.K (6) 
A + l/(r) iv + C/|*o|'(2w.-u.) = {^i-E,)v, (7) 



2M 

(^'o is taken real), and are normalized by the condition 
dY{u^ul, - Vyvl,) ^ 5^^,. (8) 



The Hamiltonian (^) does not contain the recently dis- 
cussed term originating from the presence of the "mo- 
mentum" operator of the condensate pO[ , since this term 
does not affect the elementary excitations. 

Eqs. ^ and represent a complete set of equa- 
tions for finding the wavefunctions Uu,Vi, and energies 
of the excitations. We will discuss the case of repul- 
sive (a>0) interparticle interaction in the Thomas-Fermi 
regime {fi~ nomU ^ tkoi, rtom is the maximum conden- 
sate density), where the presence of a small parameter 



(9) 



1/3n 



allows us to simplify the equations for the 



elementary excitations. First, we write Eqs. (^), and 
in terms of dimensionless eigenenergies Ei, = E^/hU 
and coordinates yi — ri/k, where li = {2^/mujfy^'^ is the 
characteristic size of the condensate in the i-th direction: 

- C^Au^ + y^u^ + {2u^ - v^)nQ = (1 + 2Ce^)u^ (10) 
-(^Kvy + y^v^ + (2v^ - u„)no = (1 - 2Ce^)v^ (11) 
-C^A*o + 2/^*0 + "0*0 = *o- (12) 

Here A^J2^i^^/^?^'^/^yl and y^ = Y.iy1- With the 
dimensionless condensate density fiQ — |^'o(r)p/nom from 
Eq.(|l|), Eqs. (|l^ and (0) are reduced to the fourth- 
order differential equations for the functions = u^±v^: 



<52/»max[eC,(C/e)'/'], 



(15) 



we can omit all terms proportional to in Eqs. (|l^), 
(|l^ ) and use the Thomas-Fermi approximation for the 
condensate wavefunction (see ^,p2|): 



*o = V"-Om(l - y^), y <1, 



(16) 



following from Eq.(^) in which the kinetic energy term 
C^A^E'o is neglected. Then, using the substitution 
f±iy) — C±(l — y^)'^^W{y), we obtain the equation 

GW + 2e^W ^0, (17) 

where the operator G is given by 

G={l-y^)A-2Y,y^i^^/oj?^/^y^. (18) 

i 

The relation between the normalization coefficients C+ 
and C_ follows from Eqs. (0), (|ll]), (|l|) and (|l8|): 



G- = eCG^ 



(19) 



The solution (|T6|) can be used in Eqs. (|Tg) and (|T^) 
from the very beginning for finding the wavefunctions 
and spectrum of elementary excitations with energies 
E^ ^ huJi. However, for the excitations with energies 
comparable to the trap frequencies this would lead to 
an incorrect result. Moreover, such a procedure makes 
Eqs. ([lo| ) and ([l^) incompatible with each other. The 
physical reason is that the wavefunctions of such exci- 
tations vary over a distance comparable with the size of 
the condensate. Hence, the kinetic energy of the conden- 
sate, omitted in the derivation of Eq.([16|), and the kinetic 
energy of the excitations are equally important. This 



is taken into account in our derivation of Eqs. (13) and 
(|r^), relying on the exact expression for ^'g. In principle, 
the exact equations (|l^) and ( p^ can be used to obtain 
a systematic expansion of the excitation wavefunctions 
and energies in the ^ parameter. 

In the case of spherical symmetry (tOi — ZJ — lu) the 
excitations are characterized by the orbital angular mo- 
mentum I and its projection m. The solution of Eq.(|T7|) 
has the form W = x''^'^P{x)Yim{9,(l}), where Yim is a 



(l-y2) -A/+ + /_ 



A*o 



A\f+-3—^AU-A f^ ° 



-A(l-z/)/_ + (l-y2)/_ 



A^-ol 
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^0 J 





^o 



~ n A*n ~ 



*0 



-3A /_ 



A*o 

*0 ; 



A*o 
*o 



= 2e% 



--2e^f- 



(13) 



(14) 



Here we have omitted the index v and written the terms 
proportional to separately. 

The low-energy excitations {E ^ or e( 1) are 
primarily localized inside the condensate spatial region. 
At characteristic distances from the condensate boundary 



spherical harmonic, x — y^, and the radial function P{x) 
is governed by a hypergeometric differential equation 



f-(¥-[)-M20) 



2 



The solution of Eq.(|20|), convergent at x = 0, is the hy- 
pergeometric function which converges at a; — > 1 only 
when reduced to a polynomial. This immediately gives 
the energy spectrum: 

Eni = hujEni = hLu{2n'^ + 2nl + 3n + /)^/^, (21) 

where n is a positive integer. The solutions of Eq.(pO[) are 
classical Jacobi polynomials P^'^^^^'"' (1 — 2y^) and, with 
the normalization conditions and (p^), we obtain 

ll/2 

/+ = 



'(l-2/2)(4n+2?+3)' 



.f-= 



£„,C(4ri + 2/ + 3)' 



^3(l-y2) 



1/2 



yipm/m(^l_2y^)YiM) (22) 
y/pO+i/2,o)(i_2y2)l^„(0,0), (23) 



where Ic = (2/i/Mu;^)-'-/^ is the size of the condensate. 
The spectrum ( |2l| ) coincides with that found by Stringari 
from the analysis of the density fluctuations in the 
hydrodynamic approach. 

In the non- symmetric case with tui ^ loi ^ ui^, the 
operator G is invariant under the inversion of any of the 



(1 



(2|m| + 1) d 
P dp 



three spatial coordinate. Therefore, the polynomials W 
determined by Eq. (^ can be labeled by the correspond- 
ing parities Vi = [±,±,±]. Another quantum number is 
the order N of the polynomial W . For even TV the func- 
tion W contains the powers of y equal to iV, — 2, 0, 
and for odd N the powers are equal to TV, — 2, ...1. 

In fact, there are only two independent parities, since 
WiVi — (—1)^. The first few eigenstates can easily be 
found. For A = 1 we have three eigenstates describing 
the condensate center of mass oscillations with the trap 
frequencies. Accordingly, there are three eigenfunctions 
W (X yi with parities Vi ~ [—] . The corresponding eigen- 
frequencies are J7 = uJe = Wi. In the case of A^ = 2 we 
again obtain three eigenstates corresponding to the con- 
densate center of mass oscillations. The eigenfunctions 
are W oc yiyj {i ^ j), the parities Vi = Vj = [— ], and 

In addition, there 



the eigenfrequencies fl — , , ~j ' ~j- 
are three eigenstates with N = 2 and parities Vi — [+] 
for all i. Those correspond to the quadrupole oscillations 
of the condensate, the center of mass being at rest. The 
eigenfunctions can be written as 

3 

W (xl+Y,b^{u}/'^^?yf■ 
i=l 

There are three different sets of coefficients bi correspond- 
ing to the three eigenfrequencies determined by the 
secular equation det[iS'] = where S is the 3x3 matrix 



The coefficients bi are determined by the system of three 
linear equations. The first one is + (^^^/w^) = 0. 

The two other equations can be any of the three linearly 
dependent equations ^ijbi = 0. 

For cylindrically symmetric traps {lui — L02 = Wp, 
C1J3 = uJz) the projection of the orbital angular momentum 
on the z axis, m, is a conserved quantity. The eigen- 
states of the excitations can be labeled by the quan- 
tum numbers A^, m, and the axial parity Vz- The 
radial parity describing the behavior of the eigenfunc- 
tions with respect to simultaneous inversion of the two 
radial coordinates is for even m, and [— ] for odd 
TO. The polynomials W can be represented in the form 
W = pl™li3„rn(p, z) exp {im(j)), where p and z are the di- 
mensionless axial and radial coordinates, and -B„m poly- 
nomials of power n ^ N — m. Each term of the polyno- 
mial Bnm has the form 2"^p"'', where Uz and rip are pos- 
itive integers. The sum + Up takes the values 0, 2, ...n 
for even n, and 1, 3...n for odd n. The integer Up is even, 
Uz being even for even n (Vz = [+]), and odd for odd n 
(Vz = [— ])• The polynomials Bnm(p,z) and eigenen- 



ergies En 



can be found from the equation 



Bnm{p, z)=0, 



(24) 



which follows directly from Eq. (pT[) . The quantity (3 — 
ujzl^p is the ratio of the axial to radial frequency. 

The number of eigenmodes fc at given vn and n is de- 
termined by the quantum number n. Generally speaking, 
we obtain k = {n + 2)/2 for even n, and k — (n -f l)/2 
for odd n. The eigenenergies can be found from the k-ih. 
order secular equation det[S'] = 0, where S is the three- 
diagonal k X k matrix 

S,j^(^~\m\\s,,-2(k-j)(k-i + \m\ + m, + S,^i,,) 



-(3^[(i^l+q)(2i~l)S^j + i(2i~l + 2q)S,+i.jl 



(25) 



and i,j — l,2...fc. The coefficient q — for even n, and 
q = 1 for odd n. 

For n = (|to| > 1) we have purely radial oscilla- 
tions, with -Bom — const and fJom = \/\n^\i^p- The 
case n = 1 corresponds to the radial oscillations, in 
combination with the axial oscillations of the center of 
mass of the condensate. Here we have Bim oc z and 
film = y^lmjwp -I- w^- In both cases the coupling between 
the radial and axial motion is absent, and the condensate 
frequencies riom and Qim are the same as those for a clas- 
sical gas in the hydrodynamic regime (see |16 !%)■ 



Si- 



l + {2-n^/u^)6,,j. 



For n > 2 the coupling between the radial and axial 
degrees of freedom becomes important, and the conden- 
sate oscillation frequencies will be different from the fre- 
quencies of a classical hydrodynamic gas. If n ~ 2 there 
are two coupled shape oscillations of the condensate with 
frequencies ftfm given by 



3 



2\m\+2 + -(3^± 



( I m|+2--/32)2+ 2/32(1^1+1) 



1/2 



(26) 



In the simplest case of to = Eq.(26) gives the fre- 



quencies of the quadrupole shape osciUations of the con- 
densate. The frequencies flfg and the frequency of 
quadrupole radial oscillations f2o2 were found in the hy- 
drodynamic approach in They were also obtained 
in |]l5| by considering the condensate evolution under a 
weak modulation of the trap frequencies, in |lq| from 



the Hamiltonian of the scaling dynamics, and in |15|| on 
the basis of variational approach. The frequencies rio2 
and ri^Q were measured in the JILA experiment |^ for 
P — \/8 and calculated numerically for this trapping ge- 
ometry in [l0|] . The frequencies were found in the 
MIT experiment for (3 = 0.08. 

For n = 3 we have two coupled shape oscillations which 
are now also coupled to the oscillations of the center of 
mass of the condensate. In this case we obtain 



3m _ 



2\m\+2 + ^p^± 



(|to|+2- 



1/2 



(27) 



Interestingly, for certain values of the projection of the 
orbital angular momentum, to, we find an accidental de- 
generacy in the spectrum of excitations. The simplest 
example concerns the frequencies ftim = u}p^\m\ + (fi 
and ^2m- As follows from Eq.(p6|), — f^'' '^'^ 
projection m satisfying the condition 



0^ - |m| + 3 



(28) 



and, accordingly, integer > 3. For /3 — \/Z we have 
TO = 0, i.e., the frequency of quadrupole shape oscil- 
lations fi^o coincides with the frequency of axial oscilla- 
tions riiQ. In the JILA trapping geometry, where /3 = -v/S, 
Eq.dl) gives \m\ = 5. 

Although the condensate frequencies do not sig- 
nificantly differ from those for a classical gas in the hy- 
drodynamic regime (for m = see p^ ), the accidental 
degeneracy determined by Eq.(p^) is characteristic only 
for the condensate. The presence of the accidental degen- 
eracy can strongly influence the picture of the condensate 
oscillations. The coupling between the degenerate modes 
is provided by the interaction terms in the Hamiltonian 
(|l|), proportional to and and omitted above in 
the derivation of the Bogolyubov-De Gennes equations. 
Therefore, driving only one of the degenerate modes, it 
is feasible to expect the appearance of oscillations rep- 
resenting a superposition of the two modes. This phe- 
nomenon ensures a clear distinction between the conden- 
sate oscillations and the oscillations of a classical hydro- 
dynamic gas and, hence, can be a signature of BEG for 
the gas in the hydrodynamic regime. 
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